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Abstract 


An efficient full-wave Green’s function technique is presented for the analysis of 
the multilayer dielectric resonator antenna (DRA) structure with an arbitrary 
number of layers. Three planar/non-planar feed configurations are investigated, 
viz. the microstrip aperture-coupled feed, the probe-coupled feed and the 
rectangular waveguide fed multilayer DRA. The modes of the multilayer DRA can 
be identified from the analysis which can be used to explain the broadband nature 
of the coupling. The layer permittivities and thicknesses are optimized for 
broadband operation of the antenna structures. The mode-excitation for the 
antenna structures is thoroughly characterized using the Green's function 
approach, providing further insight into the radiation mechanism and bandwidth 
extension for such structures. Frequency tuning of the antenna structures is also 
demonstrated by exciting the antenna at a higher order mode, maintaining the 
broadband characteristics. 

Keywords: DRA, multilayer, Green's function, antenna feed, broadband, slot- 
coupled, probe-coupled, waveguide 
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Geometry of the microstrip slot coupled N-layer DR Antenna 


(a) Side view, (b) Top view 
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Equivalent Problem for the multilayer Hemispherical DRA 


Multilayer DRA Green’s Function 
Ну 


Gy y 
Gin, =Gp + Gy 


+ Particular Solution alone represents the source radiating in an 
unbounded dielectric medium and is similar to the single layer 
DRA 


e Homogeneous Solution accounts for the dielectric 
discontinuities 


= Mode matching method is employed to analyze the Multilayer 
DRA 


1) M, excites TE to r mode 
2) M, excite both TE to r and TM to r mode 
*Three potential Green's functions are required 
1) Gy, denotes electric potential due to M, 


2) GẸ denotes electric potential due to M 
M, 4 Ld 


3) бу, denotes magnetic potential due to M, 


Potential Green’s Functions for Particular case 
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Fig. 6: Equivalent Problem for Particular Solution 
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Evaluation of unknown constants ... 


From Harrington, p 269 
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Evaluation of Unknown Constants 
in Particular Solution 


Boundary conditions in Particular 


d solutions are 


Ej, = Ер at r=Pr 


Ед –Ер =М, at r=r 


Evaluation of Unknown Constants in Particular Solution 


Orthogonal Property 
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1 Wronskian relation 
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Potential Green’s Functions for Homogeneous case 
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Evaluation of Unknown Constants in Homogeneous Solution 
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Applying the Boundary Conditions, tangential components of both E & H 
will be continuous at the dielectric interfaces at. p = НИ 
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Orthogonal Properties 
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The set of equations corresponding to о: аге 
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Homogeneous Green’s function is given by 
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Unknown constant in the homogeneous Green’s function 
can be obtained from the matrix equations given by 


т 
А [В b, b. b, bs c Draws) basa) Draw pr 


-[s 00. 0 xg 0 0 of 


т 
Али, [ба ба ба бш бы © Ges. Gora Coen. Gon] 


-[s 0 0 — 0 ug 0 0 — of 


24 


Where 


ATE, = 
Ту ту 
0 ma 

0 m 

0 ei 








m, 


0 


т, 


(ума) 





о 


т 


m, 


хут 


е 


(М-Іү2М-3) 


0 
0 


(NI -3) 


0 


My yan-2) 


Туом-э) 


0 
0 


XP yan) 


XyMy 2y- 





MM yayan-y 


My ans) 


Хт ламу 


Ty (oi) 








0 
Myan) 


0 
0 


0 


ol y (a) 


25 





0 ace 0 


ит 


путу 





0 


my 


m, 


WAN 2N-4) 


0 
0 





пут, 


ET 


0 


ШЕТЕЛДЕ 





Mayans) 


My IM yaya) 


0 


4) 


oo 


Pav) 


Myan- 


UM yiyan- 





HATTE NEA 





2) 





Plagas) 


Myan- 


HATH yo) 


UyMy any 


шт, 


NON) 


26 


мћеге 


= ЈАК), тр=—Ј(бл), т, ==, (ел), т,=Ј(Њљ), 


т, = N (Ks), m, 7, (ks), ms —-R, (s). myner- я, (аб), 


Mynx- =Å, (kuata), Manana ==, (kutsa) m, ai = =, (Вујеа), 


=F, (kytie)s My ancy = Ny (Kyte) ту o - RU (rs). 





i=1,2,3,---N 





27 


Characteristic equations of N-Layer DRA 


* Characteristic equation of N-Layer DRA for TE to r modes can be 
obtained by taking determinant of matrix ATE, and setting it to the zero. 


Пе(АТЕ,)-0 


* The Characteristics equation of TM to rmodes of N-layer HDRA can be 
obtained by taking determinant of matrix ^7M, and setting it to the zero. 


Det(ATM ,) -0 


* The Characteristics equation is independent of m. 
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* Roots of the characteristic equations corresponding to the TE to r 
and TM to r modes are complex. 


тр Z Pump — J Pump 


* Real part of the root gives the resonant frequency of N-layer HDRA 


* Real part and imaginary part of the root can be used to obtain Q 
factor of the given mode 


Onn stored ener; 
О, = 8 = ip ОШУ 
20, radiated energy per cycle 


птр 
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Mode designation 


The first subscript n in the mode designation refer to the order 
of the spherical Bessel function (Schelkunoff type ). 


The second subscript m in the mode designation refer to the 
azimuthal variation of the field. 


The third subscript p corresponds to a discrete set of roots for 
each integral value of n >=1. 
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Three layer HDRA on the microstrip substrate 
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Fabricated Hemispherical DRA and Hemispherical shells 


using Eccostock material of relative permittivity 9 and 4. 
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Stub Resonance 


3.42 GHz 





Normalized Input impedance 


be 
N 


Reactanse 


v 


285 3 38 4 45 5 55 6 
Frequency (GHz) 





55 6 





38 4 45 
Frequency (GHz) 


(а) (b) 





Normalized Equivalent Series Impedance 


^*% 


280% 





Fig. 9: Calculated, simulated and measured return loss and normalized input impedance of the microstrip 
slot-coupled three-layer DRA with centred slot excited in the 7E, тоде : 2,=9. 6, =1, 25=4, 








75mm,  n-225wm.  L,-l4Sum, W, = 1.0m, 


58mm. Ly =6.25mm, Al=0.25nm (a) Return loss. (b) Normalized input 





Omm, у, =0ит, 





т, h- 


6, =22,1=0.508 тт. TF, 





impedance. (c) Normalized equivalent series impedance 2, 








First resonance dip at 3.08 GHz is close to the source free resonant 
frequency of TE,,, mode at 3.18 GHz. 


Second resonance dip at 3.42 GHz is stub resonance. 
Calculated 10 dB impedance bandwidth is 26.80%. 
Measured 10 dB impedance bandwidth is 29.17%. 


Antenna is easier to fabricate due to the air in the second layer. 
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Frequency Tuning 
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(a) (b) 
Fig. 10: Calculated and simulated return loss and normalized input impedance of the microstrip slot- 
coupled three-layer DRA with centred slot excited in the TE, mode : r =18.5mm, 


23.5mm „L, =12.0mm , L} =4-45mm . The other parameters are the same as in Fig. 9. (a) Return loss. 


(b) Normalized input impedance. 
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First resonance dip at 4.02 GHz is stub resonance. 


Second resonance dip at 4.88 GHz is close to the source free 
resonant frequency of TE,,, mode at 4.74 GHz. 


Calculated 10 dB impedance bandwidth is 30.18%. 
The center frequency for the return loss bandwidth in Fig. 10(a) is 


at 4.44 GHz which is displaced from the center frequency of 3.36 GHz 
in Fig. 9(a) by 32.14 %. 
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Fig. 11: Calculated, simulated and measured return loss and normalized input impedance of the microstrip 


slot-coupled three-layer DRA for y-offset slot : 5 





65тт, қ-22.0тт, 1,-12.0тт. 775mm , 





Y» =4mm. The other parameters are the same as in Fig. 9. (a) Retum loss. (b) Normalized input 


impedance 
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First resonance dip at 4.56 GHz is close to the source free resonant 
frequency of TE,,, and TE,,, modes at 4.63 GHz. 


Second resonance dip at 5.54 GHz is near the loaded slot resonance 
at 5.28 GHz. 


The third resonant dip in the return loss characteristics at 6.48 GHz in 
Fig. 11(a) is close to the degenerate TE,,, and TE,,, mode resonances of the 


three-layer HDRA at 6.26 GHz, which will be examined later. 


Measured 10 dB impedance bandwidth is 43.42%. 
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Slot offset along slot width 
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Fig. 12: Calculated and simulated return loss and normalized input impedance of the microstrip slot- 


coupled three-layer DRA for x-offset slot : L, =10.5mm, Ly =3.85mm , x, =4mm . The other parameters 


are the same as in Fig. 9. (2) Return loss. (b) Normalized input impedance. 
39 


* The first resonant dip at 4.44 GHz is caused due to the excitation of 
the 7M,,, mode with a source-free resonance frequency of 4.55 GHz. 


* The second return loss dip at 4.92 GHz corresponds to the loaded 
resonance of the TE,,, mode with a corresponding source-free 


resonance at 4.93 GHz. 


* Computed 10 dB impedance bandwidth is at 21.52 96. 
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Radiation Pattern of the three-layer HDRA fed by 
centred microstrip slot 
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Fig. 14: Calculated, measured and simulated radiation patterns for centred 
slot at 3.36 GHz corresponding to Fig. 9. (a) x-z plane. (b) y-z plane 
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Radiation Pattern of the three-layer HDRA fed by 
y-offset slot. 
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Fig. 15: Calculated, measured and simulated radiation patterns for y-offset 
slot at 5.5 GHz corresponding to Fig. 11. (a) x-z plane. (b) y-z plane. 
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Maximum Bandwidth Realization of microstrip 
slot coupled three-layer HDRA 





8 





‘Maximum Bandwidth (%) Near TE, y, Mode 


* 


(a) (b) 
Fig. 16: The calculated Maximum optimized percentage bandwidth near TF,,, mode of microstrip slot 
coupled three layer HDRA as a function of ¢,, for different values of £,,. 1 =12.5mm, т, =175nm, 
қ-225тп, НИ; -10тт. х,-Отт, у,-0ит, 6,-22, һ-0.508шт, №’, =1.58т . (а) Maximum 


bandwidth versus e, (b) Source free Quality factor of TE}; mode versus £, 
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The mode excitation versus frequency for the three- 
layer hemispherical DRA excited by a centrally located 
y-directed magnetic point source 


|ТЕ/ТМ ар 


Maximum Relative Power (dB) 








Frequency (GHz) 


Fig. 17: Excitation efficiency of modes with frequency for the microstrip slot-coupled three-layer DRA 
for centred slot: £, =9, £, =1,€,,=4, 4 =12.5mm , ғ, -175иш, қ 





22.5mm , x, = 0mm , y, =0mm . 


B. Ghosh and A.B. Kakade, “Mode excitation in the microstrip slot-coupled three-layer 
hemispherical dielectric resonator antenna,” IET Microwaves, Antennas Propag., vol. 10, pp. 
1534-1540, Jul. 2016. 44 





lomogeneous Green's function of a multilayer DRA is independent of the 
modal index p, and as such the radiation efficiency of a mode also does not 
depend on p. 


»The excitation of the TE,,, modes is stronger by more than 95 dB relative 
to the other modes. 


a There is no contribution to the admittance term Y, from TE modes when 
n + mis odd. 


= There is no contribution to the admittance term ve from TM modes when 
n + mis even. 


a Thus, modes for which n + m is even can be classified as the TE type and 
for which n + m is odd can be inferred to be of the TM type. 


“ Тһе contribution of TM modes is negligible because only r- directed 
magnetic current is present. 
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The mode excitation versus frequency for the three- 
layer hemispherical DRA excited by a y-directed 
magnetic point source placed at y-offset 
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Fig. 18: Excitation efficiency of modes with frequency for the microstrip slot-coupled three-layer DRA 
for y-offset slot: 2; 79.5, 71.6, 74. я =12.5нии, б =163тт, 1 =22mm , x, = Omn , yy = 4mm . 
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+ From figure it is seen that the contribution from TE modes with n+m even 
is significant. 


“ Тһе contribution of TM modes is negligible because only r- directed 
magnetic current is present. 
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The mode excitation versus frequency for the three- 
layer hemispherical DRA excited by a y-directed 
magnetic point source placed at x-offset 
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Fig. 19: Excitation efficiency of modes with frequency for the microstrip slot-coupled three-layer DRA 
for x-offset slot: £, =9,8,3=1,8,, =4, 4 =12.5mm, n 








17.5mm , r, =22.5mm , x, = тт, у, = Отт. 
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+ From figure it is seen that the contribution from TE modes with n+m even 
is significant. 


* The contribution of TM modes with n+m odd is significant 
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Analysis of Co-axial probe coupled three-layer HDRA. 
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(b) 
Fig: 22: Geometry of the probe coupled multilayer DR antenna. (a) Top view. ( b) Side 
view 
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Probe fed HDRA problem is dual of slot fed HDRA. 


The centred slot would couple to TE,, modes ofthe multilayer 
HDRA. 


The centred slot fed HDRA produces broad side radiation pattern. 


The centred probe would couple to 7m,,, modes of the multilayer 
HDRA. 


The centred probe fed HDRA produces monopole like radiation 
pattern due to excitation of TM;,, modes of the DRA. 
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Fig. 23 : Calculated, simulated and measured return loss and normalized 
input impedance of the probe-coupled three-layer DRA: ¢,, =9, €,. =1 , 
2454, пе125тт, т, =140тт, r, =20.5mm, 1=6.7mm, 4=0.63mm, b= 2.0mm, 


c-0mm , (a) Return loss. (b) Normalized input impedance. Es 


* The first resonant dip at 4.44 GHz is caused due to the excitation of 
the 7M,,, mode with a source-free resonance frequency of 4.48 GHz. 


* The second return loss dip at 5.64 GHz corresponds to the loaded 
probe resonance at 5.44 GHz. 


* Measured 10 dB impedance bandwidth is at 43.796. 
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Radiation Pattern 
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Fig. 24 : Calculated, measured and simulated radiation patterns for the 
centred probe fed DRA at 5.14 GHz corresponding to Fig. 23. (а) х-2 plane, 


(b) y-z plane 


Multilayer HDRA excited in 7/,., mode produces monopole 


like radiation pattern 
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Maximum Bandwidth Realization of Probe 
coupled three-layer HDRA 
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Fig. 26: The calculated Maximum optimized percentage bandwidth near 7M,,, mode of 
probe coupled three-layer HDRA as a function of £s, for different values of 
215 =12.5тт, n =14.0mm, n =20.5mm, c=Omm. (a) Maximum bandwidth versus 


6, 


£,,. (b) Quality factor of TM,,, mode versus ¢,;. 
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Wide band operation 
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Fig. 27 : Calculated return loss of the probe-coupled three-layer HDRA 
for wideband operation: ¢,, =9, 7, =12.5mm , r, - 14.0mm, r,=20.5mm, a =0.63mm, 
Ь-20тт,с-0тт (а) 6, = 75, 1-5.8тт, ВМ/ - 63.68%, 
(b)5;-4 54-3 I= 6mm, BW = 63.22%. 
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Analysis of the Rectangular Waveguide Slot Coupled 
Multilayer HDRA 
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Fig. 28; Multilayer hemispherical DRA slot-coupled to shorted end of rectangular 


waveguide. (a) Side view. (b) Front view 











Multilayer DR 


Waveguide Antenna 








(b) 


58 





A.B. Kakade and B. Ghosh, “Analysis of the rectangular waveguide slot coupled multilayer 
hemispherical dielectric resonator antenna,” JET Microwaves, Antennas Propag., vol. 6, pp. 338- 
347, Feb. 2012. 59 


Summary of the problem 


The scattered magnetic field inside the waveguide due to an 
aperture in the transverse plane is derived. 


The ground plane thickness is taken into account in the analysis 
by modelling the slot thickness as a rectangular cavity. 


The Green’s function of a rectangular cavity due to magnetic 
point source is derived. 


In the MoM analytical integrations are performed inside the 
rectangular waveguide and the slot cavity, whereas DRA 
Green’s function is integrated numerically using the PTDI 
technique. 
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Fig. 29: Calculated, simulated and measured coupling characteristics of the three-layer 


DRA excited by the empty WR-90 waveguide: ¢,,=1, ¢=4. 6; =9, 2Z,=15.0mm, 
27, 


=Omm (a) Return loss. (b) Normalized input impedance. 





lOmm. n-l3.5mm.r,-l6.5mm. тһ = 20.0тт. 21=1.3 тт. х, = Отт. 
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e It has been found that wide bandwidth of the rectangular 
waveguide slot coupled to three-layer HDRA is obtained when 
the permittivity of the innermost layer is unity and the dielectric 
filling the waveguide is air. 


e The 10 dB impedance bandwidth of ~ 12.9% is achieved with an 
empty rectangular waveguide slot coupled to three-layer HDRA. 
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Fig 29 A: Calculated, measured and simulated radiation patterns at 9.86 GHz corresponding to 
Fig. 29. (a) x-z plane. (b) y-z plane. 
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Fig, 29 B: Calculated and simulated coupling characteristics of the three-layer HDRA excited 
by the empty WR-90 waveguide with slot-offset along x-axis: &; —1. $,—4. &,—-9. 
21,=14.1шш. 20, =1.0 пип, 7-135mm, 4=163 mm, 14=19.7 mm, 
3. 0 mm (a) Return loss. (b) Normalized input impedance. 
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«The first weak resonant dip at 10.04 GHz is caused due to the excitation of 
the TE,, mode with a source-free resonance frequency of 10.06 GHz. 


* The second return loss dip at 10.74 GHz corresponds to the loaded 
slot resonance at 10.64 GHz. 


* Calculated 10 dB impedance bandwidth is at 5.42%. 
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Thank 
you 


